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Abstract 

Individual generators in multi-machine networks are 
represented by closed-loop configurations of linear time- 
invariant modeling and nonlinear feedback, subject to 
modeling uncertainties and power disturbances, based on 
which control laws and implementation algorithms of static 
VAR compensators (SVC) are established to robustly 
stabilize individual generators by combining the gain- 
scheduling technique to the small-gain theorem in the sense 

of finite-gain L 2 -stability. The gain-scheduling SVC 

manipulation approach can be implemented in multiple 
decentralized algorithms that involve only locally-measured 
signals, and are robust against modeling uncertainties and 
power disturbances. The suggested SVC control laws can 
also accommodate the steady-state performance 
specification in term of generator power angle, besides 
stabilization. As a special case of the suggested SVC control 
laws, examples are included to illustrate efficacy and 
simplicity of an averagely gain-scheduling SVC control law. 
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Introduction 

Stabilization of individual synchronous generators in 
multi-machine power networks is inevitable and of 
engineering significance much higher than most 
techniques for stabilizing all generators [Alberto & 
Bretas, Cong et al., Gupta et al., Yang et al.] 
simultaneously, which usually involve complicated 
control laws and expensive communication devices, 
and high maintenance cost, though some of them 
are implementable in decentralized ways [Aboul-Ela 
et al., Damm et al., Fan & Feliachi, Lu et al., Wang et 
al., J. Zhou (2010a)]. Stabilization problems in 
synchronous generators are also conventionally 
formulated as power swing reduction, to which 



numerous papers are devoted in miscellaneous 
stability meanings, and latest results are reported in 
[Chiang & Chu, Haque (2005), Liu et al., Salam et al., J. 
Zhou & Ohsawa, J. Zhou (2010b)]. 

Robustly stabilizing individual generators in a multi- 
machine network in the finite-gain L 2 -stability via 

flexible ac transmission systems (FACTS) [Haque 
(2004), Lerch et al.] is considered in the study. Control 
laws of static VAR compensators (SVC) are developed 
by exploiting the gain-scheduling mechanism [H. K. 
Khalil(2002)] and the small-gain theorem [Vidyasagar]. 
The gain-scheduling SVC approach possesses many 
advantages, for example, it can be implemented in 
multiple decentralized fashions that involve only 
locally measured feedback signal (and thus, no time- 
delay concerns in power networks need to be taken 
into account); and it can cope with power disturbances 
and modeling uncertainties such as those in damping 
windings, and results in robust stabilization for us. An 
implementation algorithm for the suggested SVC 
control law is also contrived, which can also 
accommodate the steady-state power angle 
specification, besides stabilization. However, our 
numeric simulations show that the suggested SVC 
control laws may be conservative, compared to those 
of [J. Zhou (2010b)]. 

By the literature, neither damping windings nor 
disturbances can be dealt with accurately and 
effectively by the equal area criterion (EAC) [Haque 
(2004, 2005), E. Z. Zhou]. EAC is a graphical tool, and 
requires expertise operators with expertise. Lyapunov 
theory is also frequently employed for stabilizing 
generators, which usually yields sophisticated 
nonlinear algorithms with high conservatism [Yang, J. 
Zhou (2010a)]. Recently, pole assignment and least 
quadratic regulation are adopted for SVC control law 
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design [J. Zhou (2010b)], whose engineering efficacy 
depends closely on modeling, measuring devices and 
actuators. 

The second section collects facts on swing equations of 
individual generators with SVC's connected, their 
state-space re-expression and feedback configuration 
interpretation; and stabilization problem is formulated 
as well in the section. The third section introduces 
concepts and useful stability criteria. The next section 
talks about stabilization control laws of SVC's by 
combining the gain-scheduling technique and the 
small-gain stability theorem. The third and fourth 
sections are the main context. Illustrative examples are 
sketched in the fifth section with an implementation 
algorithm of the average gain-scheduling SVC control 
law, while the final section states our conclusions. 

Preliminaries and Problem Formulation 

Individual Generator Systems and Equivalent Swing 
Equations 
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FIG. 1 IGS CONFIGURATION IN A MULTI-MACHINE 
NETWORK 

An individual generator system (IGS) is isolated from 
the multi-machine network as in Figure 1, where 

IGS represents the remaining multi-machine portion 
connected with the IGS via the transmission line (TL). 
B g , B m and B c , respectively, are the internal bus, the 
terminal bus to which a shunt SVC is connected, and 
the local load bus to which TL is linked. X is the 

reactance between B m and B g , while x l is that 
between B m and B c . SVC is parameterized as 
variable reactance X svc as in [Haque (2004, 2005)]. 
P m is the mechanical power from turbine and 
governor, and P is the electrical power viewed at 
B g . CL represents the control law that manipulates 
the SVC. 

Based on [Haque (2005)] and the re-modeling 
technique of [J. Zhou (2010b)], an equivalent swing 
equation depicting the dynamics of the generator rotor 



with a shunt SVC as in Figure 1 is 
8 = CO 



CO = 



J 



co- 



b sin 8 P + P. 



(1) 



xM 



where b = E c E g and x = X + X, 



M 

X svc X g X, with 
E c and E are the virtual voltages of B c and B g , 
respectively. Voltage deviations from E c and E g are 
reflected by perturbation in P m . Complicated but 
weak dynamics of the SVC are viewed as uncertainties 
and added to P A . It is hard to handle P A exactly. It is 
assumed that P A is a bounded disturbance. 

Also in (1), 8 and CO are the generator rotor angular 
displacement and velocity, respectively, with respect 
to a reference axis on B c rotating at the system 

angular frequency CO^ [Kimbark]; J : combined inertia 

moment of the generator rotor (kg-m 2 ); andM = Jco s ; 

D : mechanical and electrical damping due to 
amortisseur windings [Anderson & Fouad, 
Machowski et al.]. 

State-Space Expression of IGS and Closed-Loop 
Configuration 

In what follows, let 8 be the output of the IGS, that 
is, y = 8 . In the steady state, 8 is called the power 
angle of the generator, which is related to the power 
factor that reflects the power generation efficiency of 
the generator. Based on (1), the IGS is described by the 
state-space representation 



g = Ag + B(P m +P A ) + 

y = cg 

where g '.= [8, cof and 
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FIG. 2 FEEDBACK CONFIGURATION INTERPRETATION FOR 
THE IGS 
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Examining (2) carefully reveals readily that it is 
essentially nonlinear, and possesses a feedback 
configuration interpretation of Figure 2. Now we 
briefly denote 



I GO) = C(sl 2 - A) l B = [Ms(s + D/J)Y 
[H(x, y) = x l bsin y 



(3) 



It is self-evident that G(s) is the transfer function of 
the LTI subsystem G : e l— > y and H(x,y) is a 
nonlinear function reflecting the mapping H : y I— > u . 

Finite-Gain Li-Stabilization Problem Formulation 

Denote the Hilbert functions linear space 

L 2 =: \x(t) g i?"|(JJ|x(r)| 2 (ir) 1/2 < oo . Finite-gain L 2 - 

stability definition and the small-gain theorem can be 
found in [Khalil, Vidyasagar]. Under the configuration 
of Figure 2, the small-gain theorem [Khalil, p. 218] 
provides us with a sufficient condition for the closed- 
loop system to be finite-gain L 2 -stable as follows. 



IIL2/L2 



W\\ <i 

II Wliili 



(4) 



where 



is the L 2 -induced norm of (•) 



Wliili ~ ~' 2 
Stabilization in finite-gain L 2 -stability of the IGS is 

formulated as follows: determine SVC reactance 
manipulating laws such that (4) is satisfied. 

Gain-Scheduling Small-Gain Theorem 
Interpretation 

Loop Transformation 

One might suggest that we use the H x norm of 

in (4). We have 



G(s) in place of G| 

IIg" 



Li/L, 

= \\g\\ 



\\GU<P)\\ 2 



Wliili • su P<p^(-« 

= ™p <pe{ _^ } \G(j(p)\ 

However, G(s) possesses a pole at the origin, and 

thus G I = co . This hampers us from exploiting the 

small-gain inequality (4) directly in determining 
control laws to SVC. 

To surmount |G||^ = co , a loop transformation K is 

introduced to the configuration of Figure 2 as in 
Figure 3. The input/output relationships in Figure 2 
and 3 are the same. It is clear that 
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H(x, -,K) : y i-> u 
FIG. 3 LOOP TRANSFORMATION IN THE IGS CONFIGURATION 

*, G(s) 1 

G(s, K) -. = 

l + KG(s) Ms(s + D/J) + K 

u=H(x,y)-Ky = (- 



Mny -x>-K)y 



y 



(6>(y)x i -K)y.= H(x,y,K)y 



(5) 



where 9(y) = bsm yl y . The LTI subsystem 
G(s, K) is stable as long as K > . 

Under the configuration of Figure 3, the small-gain 
inequality is 



G(;K) 



H(X,;K) 



<1 



(6) 



Apparently, 



G(;K) 



is equal to the //^norm of 



G(s,K) , namely G(-,JQ / which can be given 

I 1 1 oo 

numerically via the bi-section algorithm [K. Zhou & 
Doyle]. A closed-form expression for G(-,iT) is 

available in Appendix A. However, 



w(x,-,^) cannot be given explicitly due to its 

II ' \\L2IL2 

obvious nonlinear feature. 

Gain-Scheduling According to Phase Plane 
Partitioning 
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FIG. 4 PHASE PLANE PARTITIONING 
To reduce conservatism of the small-gain inequality 



(6), 



G(;K) 



and lff(x,-,2n should be 

L 2 /L, II '' lli,/L, 

estimated as accurately as possible. To this end, the 
gain-scheduling technique is employed [Khalil]. We 
firstly partition the phase plane of (2) as in Figure 4, in 
which the solid lines are the left (included) boundaries 
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for Dj and the dashed lines are their right (not included) 
boundaries; namely D ; =: { (8, CO) G R 2 : 8 j < 8 < 8 M } . 
The sub-regions are sufficiently sufficient so that 
U, D j contains all phase portraits of (2) that interest 

us, while each Z). is narrow enough (or d is 
sufficiently small) so that 8 = 8 { is reasonable for 
every 8 G D. . 

According to the phase plane partitioning, it is 
assumed that for all 8 = y G D i , 

G(;K) = G(;K i ) 

H(x,;K)y = [(0(8 i )xr l -K i ) + V i ]y 
= [H(x l ,8 i ,K l ) + V i ]y 

where K t > is a loop transformation introduced for 
8 = y G £). ; namely the loop transformation of Figure 
3 is replaced with a gain-scheduling one 
{•••,K_ 1 ,K ,K 1 ,---}—.K defined according to 

{•• •, D l ,D Q ,D l , • • •}=: D . K i is used as a design 
parameter over D { . Correspondingly, x { is the control 
action of the SVC for 8 = y G £). ; namely, the control 
law of the SVC is also gain-schedulingly manipulated, 
denoted by {•■•,X_ 1 ,X ,X 1 ,---} =: X . In addition, V ; 

means the modelling mismatch when H(x,-,K) is 
approximated with HijK^S^K^) ; clearly, V ( . 
depends on x t and 8 e D { but independent of K i . In 
the following, all V. are simply viewed as bounded 
modelling uncertainties and collectively denoted by 
V :={•••, V_ 1 ,V ,V 1 ,---};namely 



Pm+PA +B 





lim^ |V ; .| = 
sup |V ( .| < c 



(7) 



for all i . In (7) d is the sub-region width defined in 
Figure 4, and < c < co is a (small) constant. 

With the gain-scheduling approach, Figure 3 is re- 
drawn as a gain-scheduling configuration of Figure 5. 
By the notations of Figure 5, the small-gain inequality 
(6) is modified as 

P(-' K i L2IL 2^M\H(X,-,K) + V\\ L2/L2 <l (8) 
vv 




CM 



«(X,-,K) + V:y~ti 



FIG. 5 GAIN-SCHEDULING LOOP TRANSFORMATION 
CONFIGURATION IN THE IGS 

To exploit the gain-scheduling framework, we must 
answer: i ) areG(-,/0and//(X,-,Z<0 + V meaningful as 

mappings on L 2 ? ii) can we get upper bounds for 

llGO.AOL,, and \\H(X,;K) + V\\ r , , explicitly? 

I II Ln I Ln i II Ajj I L~, 

These issues will be addressed in the subsequent 
subsections. 

Well-Definedness of H(X,-, K) + V and Upper Bound 
for\\H(X,,K) + V\\ L2/L2 

To understand the following deductions, 7, is utilized 
to denote the union of all the time intervals when 
y G D i . By the definition of l t , we have 

lU=[0,co) 

Note that each H(X/ , 8 t , K ( ) + V ; a scalar constant, is 
followed by the definition that 

H(X,,K) + V h/u 

= sup\[H(x,y,K)+v]y\\ J\\y\\ 

= sup [\o\[H(X, y(r),K)+VW\ 2 dr\ m M 

= sup [J u , \[H(X, y(T),K)+V]y(TfdT] V2 M 

Mli^ ' 

= sup E I, \[H(X, y(T),K) + VjyirfdT]" 2 M 
= sup IX, I, \H(x. , 8. , K t ) + V, ' yirydrr M 
= sup|I, ^ f ,4,^.)+V. ? J ( kr)| 2 rfr] 1/2 /N| 



< sup [maxtsuiJ^^.^O + vJ }I, J,Jy(r)| 2 dr]" 2 

Wllvll <1 ' WX7. I ' 



2 / V 



v||>|| <i 



= sup [maxJsuplflX^., ^.iQ+vJ }l u , |y(r)| 2 rfr]" 2 /Ml 



(9) 



= max{sup|ff(^,<y j ,Je;)+vJ} sup [Jo bWfrfr] 1 ' 2 /II 
^axJsupjffC^.J^+vJ} 
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where j 7 Q(r)Jr denotes the integral of the 
integrand (-)(r) over l i . Noted by (7) that 

max{supH(x i ,S i .K i ) + V i } 

' VV; 

< max{ H{x i , 8 i .K i ) + sup|V,. |} 

' VV, 

= max{ H(x i ,S i .K i ) + V i + suplVJ} 



V£eD, 



< max{ H(x i , 8 i .K t ) } + c < oo 

i 

This, together with (9), reports that when we view 
H(X,;K) + V as a mapping from L 2 to L 2 , then it 
is bounded; in view of this fact, it is noted that 
H(X,-, K) + V is well-defined. 

Well-Definedness of G(-, K) and Upper Bound for 

To verify our statements, we define the following 
time-domain projection operators. 



1, f E /,. 

0, f £ /, 



By the definition, the following facts can be claimed. 

XJ Pi =1 , Pfi = P t , and P.P k = whenever 
i * k ; 

• P ( y stands for the pieces of y(= S) that 
belong to Z) , while P t e reflects those pieces of 
e defined on the time interval /, ; 

• G(; K)(P iZ ) = G(; )(P iZ ) , where G(; K t ) 
stands for the LTI system whose transfer 
function is G(s,K t ) as explicated in (5). 

From Figure 5, the time-domain relationship between 
e and y is 

y(t) = G(-,K)(e)(t) (10) 

In fact, G(-, K) represents a set of linear time-invariant 
differential equations that switch from one to another 
at isolated time points that form a set of measure zero; 
namely, G(-,K) can be interpreted as a piecewise LTI 
system, or equivalently it can be said that G(-, K) is 
linear almost everywhere in time (a.e. t e [0, oo) ). The 
differential equation theory [O'regan] reveals that 
such switchings have no affect on the dynamics of 



G(-, K) on a whole. This, together with the properties 
of the projection operators and (10), implies that 

y(t) = G(-,K)(e)(t) = G(;K)(L k P k e){t) 
= T k G(-,K)(P k e)(t), a.e.t c[0,cc) 

= H k G(;K)(P k e)(t), a.e.t c[0,cc) 

Now it is observed that 

\\G(;K)\\ = sup \\G(;K)(e)\\ L /\\e\\ L 

II IILWLj ^11 IIL 2 II IIL, 



f 00 I 

- sup [f \Z k G(;K)(P k e)(T) 

wll^ll ^ 1 *V I 



drf 2 l\\e\\ 



= sup [^f^GO^^^O'^* ( n ) 

wll.ll 7 77" J-=0 



vii«|, <i Ik 

11 lli 2 



■T k G(jcp,K k )E k Ucp)d(pf 2 /\\e\\ L2 

= sup [^-f \G(j(p,K)' E(j(p) 2 d<p] 112 /\\e\\ 

\/\\e\\ L2 <i In J - co1 Ll 

where (•) means complex conjugate transposition, 
and 



G(j(p,K) =: 



GU<P,K_ X ) 
G{j(p,K x ) 



eC 



N 



E(j(p) =: 



E (j<P) 



i A' 



where iV is the number of the sub-regions in the phase 
plane partitioning D . In (11), the Parseval theorem is 

used for the pair G(-,K k )(P i e) and 

G(j(p,K k )E k (j(p) . Here, E k (jq>) is the Fourier 

transform of P k e . E k (j(p) is well-defined since 

P k e e L 2 for e e L 2 . 

Applying the Cauchy-Schwarz inequality to the last 
equation of (11), we continue 
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G(;K) 



L2IL2 



< supf-'-f (G(j<p,KyG(j<p,K)) 
•(£(j»*£(j»)^] 1/2 /H L2 

= supt-^fd, G(j<p,K)*GU<p,K)) 
■(Z* Eijcp)* Eimdcpf 1 l\\e\\ Li 

< sup [—J (Z ( 



Ve L2 <1 2^" ' 



GO, AT) 



) 



(Z* E(j<p) ' E(j<p))d<p] 112 l\\e\ 



(12) 



2 ) 1/2 . In (12), we 



V II II oo 

• sup r z t £(j»*£(j»w /2 /|H|, 

ve i2 <i 2n J -°° ^ 

= : II G (-^)IL 

• supt-J-fa* £ t OW)Z, E^dtp] 112 l\e\ 

where \\G(-,K)\l =: (L k p(;K k ) 
have used the following facts: 

uJu E kV ( py E iV ( p) d p 

o (P,e)(r) (/>e)(r)dr = 
whenever i ^ k ; and 

[^F (Z k E k U<p))I. i E i U<P)d(p\ 

= [f(Z tJ P^(r))*Z ; 4 e (r)Jr] 1/2 

= [£ e (r)^(rWr] 1/2 = ||4 

The facts follow from the Parseval theorem and the 
properties of the projection operators. 

A 

Since G(-,K k ) < go for each k and we only consider 

the case with finite many sub-regions, therefore we 
have ||G(-,iOL 

< co . It is revealed that when 

G(-, K) is a bounded mapping from L 2 to L 2 ; or 
equivalently, that G(-, K) is well-defined. The well- 
definedness of G(;K) and that of H(X,-,K) + V 
validate our arguments about the small-gain 
inequality (8). 



1/2 



Interpreting the Small-Gain Stability Criterion with 
(9) and (12) 

With the upper bounds for \\H(X,;K) + V and 

II WL^/L^ 

namely (9) and (12), It can be 



\\G(;K)\\ 



concluded from the small-gain inequality (8) that for a 
prescribed gain-scheduling loop transformation K , if 
we can design a gain-scheduling control lawX such 
that 

llGO,^)! ■max{sup\H(x i ,S i ,K i ) + V i \}<l (13) 

then the closed-loop IGS is finite-gain L 2 stable. 

The small-gain inequality (13) is a sufficient condition. 
Examining (13) carefully, one will realize that 
existence and performance of the gain-scheduling 
control law X for maneuvering the SVC are closely 
related to the gain-scheduling loop transformation K . 
There are numerous ways to exploit the small-gain 
inequality (13) in determining appropriate gain- 
scheduling control laws in form of 
X = {••-,X_ l ,X ,X l ,---} for manipulating the SVC 
reactance. 

Stabilizing the IGS with Gain-Scheduling SVC 

In this section, we first discuss possible gain- 
scheduling control laws of the SVC when only 
stabilization is concerned. Then it is shown that the 
gain-scheduling SVC control law can also be 
employed for power angle specification, which is one 
of the most important steady-state performances for 
synchronous generators. Finally, it is verified that 
stabilization and steady-state performance 
specification in the IGS can be implemented 
simultaneously with a single gain-scheduling SVC. 

Stabilization under Average Gain-Scheduling SVC 
Algorithm 

It is clear that the small-gain inequality (13) holds as 
long as 



G(;K)\\- H(x i ,S i ,K i ) + V i 



<1 



(14) 



for each i and for all V ; ; that is, each x t is 
determined with respect to the gain-scheduling loop 
transformation K = {••-, K _ x ,K ,K l ---} such that (14) 
is true. More explicitly, we write 

6(S i )x i - l -K i +V i <1/\\G(;K)\\ 
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for each i and for all V, . Clearly, the above inequality 
can be equivalently interpreted as 

'^ i )* i - 1 -^+v i <i/||G(-,A:)|| ao , 

if 0(S i )xr 1 >K i -V i 
-diS^x- 1 +K i -V i <l/\\G(;K)\l, 
if 0(S i )x i - 1 <K i -V i 

for each i and for each fixed V ; . . More simply we 
write 



0{S i )xr l <K 

if 0(S t )x. 
diS^xr^K 
if 0(S t )x, 



-v, + i/||g(-,jo|| 

>K,-V, 
-V, -l/||G(-,Ar)|| 

<K; -V, 



(15) 



for each i and for each fixed V ( . It is not hard to see 
that for each i and any given K i > and each fixed 
V,- , if we can choose x i to satisfy both of the 
inequalities in (16), 

K i -W i -l/p{;K)\l<0{S i )x- 1 
< K i -V i+ l/\\G(;K)\\ 



(16) 



then at least one of the inequalities in (15) is true. 
Furthermore, it is noted that 



K i -V i +l/\\G(;K)\\ 
= 2/\\G(;K)\\ >0 



■(^-^-i/IIgo.js:)! ) 



It is observed that such x i always exists as long as 
0(S t )*O so that (16) is satisfied. 0{8 i )*Q can be 
guaranteed by the selection of the left boundary of D l 
so that sin 8 i I S^O . 

Remark 1 The gain-scheduling loop transformation 
K is utilized to derive the small-gain inequality (13), 
which yields the gain-scheduling SVC control law (16). 
K should be implemented nowhere in the closed-loop 
IGS; and the control law X satisfying (16) is only 
embodied. The gain-scheduling loop transformation 
K plays its role indirectly in the closed-loop IGS via 
the control law X . A different gain-scheduling loop 
transformation brings different effect into the closed- 
loop IGS. There are infinitely many possibilities to set 
Xj such that (16) is satisfied. 

Remark 2 A simple way to parameterize the SVC 
reactance is to choose the average evaluation of 



K t -V f +l/\\G(-,K)\l andtf. _v f +l/\\G(;K)\\ oo , namely 



G(;K) 



or 



In this averagely gain-scheduling SVC algorithm, 
there is no need to compute any H m norm for the 
involved LTI models. 

One can take K t > sufficiently large, and at the 

same time let |V ; | sufficiently small by partitioning 

the phase plane in the sense of d — > , see Figure 4. 
Then the averagely gain-scheduling SVC algorithm 
reasonably and approximately reduces to 



X r 1 =K i /0(S i ) 



(17) 



Our approximating x t = (K l . — V ; )0(S t ■) can be 

validated by (17) as follows. Let us return to Figure 5 
and notice that the feedback signal caused by the 
modeling mismatch V can simply be attributed to P A 
by separating V and H(X,-,K) (that is, Vy is 
treated as part of the power disturbance P A ), and then 
the previous arguments can be repeated similarly 
without taking V into account. In this way, we once 
again led to (17) eventually. In view of this, the 
stabilization control law (17) is robust against 
modeling uncertainties and power disturbance. This is 
highly expected in generator stabilization engineering. 

Transient and Steady-State Performances Specification 

At first, we explicate how the loop transformation 
parameter K (or equivalently the gain-scheduling 
control law (17)) is used to affect transient 
characteristics of the closed-loop IGS. 

According to Figure 3, the transfer function from 
P m + P A to y when the SVC reactance is set at the 
averagely gain-scheduling control law (17) is given by 



G c {s,x t ) = 

i + G(s, a:, 

G(s) = G(s) 

\ + G(s)e(d i )x^ \ + G{s)K i 

1/M 

" s 2 +(D IJ)s + KJM 

where the fact is employed that the loop 
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transformation K t can be removed simply due to the 

configuration. Clearly, the poles of G c (s,x j ) are 
given by 

\ 2 =-D/(2J)±^D 2 /(2J) 2 -K i /M (18) 

Remark 3 Clearly, for any K t > 0, Re(\ 2 ) < . It is 
also observed that the bigger K i > is, the more 

quickly the SVC control will be affected. Bearing this 
in mind, we can prescribe the gain-scheduling loop 
transformation K = {••-, K x ,K a ,K x ---} as follows; if 
the power angle 8 deviates far from the equilibrium 
point, then bigger K t > is taken so that the SVC 
drives the 8 -dynamics more quickly back to the 
equilibrium; otherwise, smaller K i > is taken to 
avoid overshooting. 

Secondly, steady-state performance specification is 
taken into consideration by means of the loop 
transformation parameter K (or equivalently the 
gain-scheduling control law (17)). 

Power angle of the IGS in the steady state, 
y ss = lim y(t) , is closely related to power generation 

t— >CO 

efficiency of the concerned generator. Hence, how to 
adjust the SVC reactance for the closed-loop IGS to 
attain a desired power angle performance, say S S5 , is 

significant. In this section, it is explicated that such 
SVC control laws do exist. 

Again from Figure 3, the transfer function from 
P m + P A to y when the SVC is set at x i is 



G c (s,x,) = 



1/M 



s 2 +(D/J)s + 0(8 i )/(Mx i ) 



(19) 



The steady-state response of G c (s,x i ) when P m + P A 
is a step signal in form of U ■ l(t) is 

_ hms U_ l/M _Ux L _ 

y * s ~ ™ S 7 s 2 +{DIJ)s + e(8 i )l(Mx l ) ~ 0(3) 

Remark 4 From the above steady-state relationship, it 
is seen that one can specify the steady-state 
performance of y to the desired position y ss = 8 SS by 

setting the SVC reactance x i as follows. 



Congruence between Stabilization and Steady-State 
Specification in the Same Average Gain-Scheduling 
SVC Algorithm 

In the preceding subsections, stabilization, transient 
response features and steady-state performance 
specification for the closed-loop IGS are tackled which 
are governed by (17), (18) and (20), respectively. Since 
these three objectives are expected to be attained by a 
same and single SVC, such a question is inevitable: 
whether a gain-scheduling loop transformation K 
and a gain-scheduling control law X exists such that 
(17), (18) and (20) are satisfied simultaneously. Here 
only the parametrization congruence issue between 
(17) and (20) is in consideration. 

When the IGS runs into steady state, it is assumed that 
the SVC control law is determined under a loop 
transformation K t = K Q over a sub-region Z> = D () , 

in which the desired power angle 8 SS is located. 

Hence, the congruence issue can be formulated as 
follows: to achieve the steady-state performance 
specification x i = X Q of (20), that is, 



(21) 



while a loop transformation K Q > can be fixed such 
that 



K ' ' ~ V " ~ 1 < 0(8 O ) V < " V ° + 1 



L v 

\\G(;K)\\ 



(22) 



(20) 



where the assumption is that K is sufficiently larger 
than V . 

Now we observe by (21) that 

0(S o )x o - 1 = 0(S o )&-0(S o )r l = ^- > 

Therefore, if let K =U I S ss , the steady-state 

performance specification control law (21) satisfies (22) 
immediately. This leads us to the desired congruence 
between (17) and (20). 

The congruence between stabilization and steady-state 
performance specification indicates that the steady- 
state performance specification control law of the SVC 
automatically guarantees finite-gain L 2 -stability in 
the closed-loop IGS, as long as control laws over other 
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sub-regions run according to (17) that ensures the 
small-gain theorem condition (16). 

Numeric Illustrations 

Numeric illustrations show that the averagely gain- 
scheduling control laws of SVC's given in (17) is a 
simple and effective technique to stabilize individual 
generators, while steady-state features of the closed- 
loop IGS's can be specified. Table 1 lists the 
coefficients of an example IGS. 

TABLE 1 COEFFICIENTS OF THE EXAMPLE INDIVIDUAL GENERATOR 



b : maximum power capacity 


950(kw) 


D : amortisseur damping constant 


95 


J : combined inertia moment 


550(kg.m2) 


C0 S : system angular frequency 


2 n x 60(rad 
/sec) 


P : mechanical power 


475(kw) 


%g : reactance between B m and B g 


0.6(0) 


X[ : reactance between B and B 


0.4(D) 



In the numeric simulation, all time response curves are 
plotted via the second-order Runga-Kutta algorithm 
[Vidyasagar] with [S(0),co{0)] T randomly given within 
£(0)e[-l,l] and cj(O) e[-l,l] , and the step length 
O.Ol(sec). The solid lines represent the time behaviours 
of the generator rotor with SVC installed, while the 
dash-dotted lines stand for those without SVC. The 
disturbance P A is white noise with \p\ < 20 (kw). 

SVC Actions Based on the Averagely Gain-Scheduling 
Control Law 

Only the control law suggested in Remark 2 is in 
consideration. We partition the phase plane into 
equitable sub-regions as follows: the sub-region 

D = {(S,a))&R 2 : S ss -AS<S< S ss + AS} 

where AS means a deviation from 8 . In the 
following, we take AS = 8 SS x 10% . 

All other sub-regions locating to the left or the right of 
D Q , namely Z). , are 

D. = { (5, co)eR 2 : S ss + (2i - 1) AS <S< S ss + (2i + V)AS} 

where 1 ~ ' ' _ ' _ ' ' ' " . According to the phase 
plane partitioned, a gain-scheduling loop 
transformation K = {••■, K x ,K Q ,K X ---} is prescribed. 



Taking the steady-state specification lime? — S ss into 

consideration, and for numeric simplicity, we will take 
V 



K Q = — >0,-K_ 2 =K^ -K i -K 2 



Consequently, the averagely gain-scheduling 
algorithm of the SVC is 



S„bsin(S„-AS) 
x n =— , VJgZ)„ 



U(S SS -AS) 
bsm(S+(2i-l)AS) 



K~ (S ss +(2i-V)AS) 
Simulation Results and Observation 



V£eD,(/*0) 




a 5 ia is a is- sii is m is sii 

tm t liicj 

FIG. 6 TIME RESPONSE CURVES WHEN S = 30° (OR tt/6 (rad)) 

In Figure 6, the steady-state power angle is specified at 
S = 30° or 0.5236 (rad), which is also the intrinsic 

steady-state power angle of the individual generator 
(that is, when the SVC is not installed). By this 
structural feature, it is not strange to see that the SVC 
is actually switched off in the steady state. 




10 1^ Z9 



FIG. 7 TIME RESPONSES CURVES WHEN § = 72° (OR 1.2566 
(rad)) 
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In Figure 7, the steady-state power angle is specified at 
<J 5 =72 or 1.2566 (rad). Compared to the case of 

Figure 6, since the steady-state power angle is 
different from the intrinsic power angle, the SVC must 
provide manipulation affect even in the steady state. 

The numeric results of Figure 8 clearly reveals that for 
some initial conditions, the generator rotor may run 
simply out of synchronism without enacting the SVC; 
in other words, the suggested gain-scheduling SVC 
control law does stabilize the IGS system, while the 
steady-state power angle is driven to the desired 
position S a = 72° or 1.2566 (rad). 

Conclusions 

Robustly stabilizing individual generators in multi- 
machine networks in the finite-gain L2-stability sense 
via output-feedback-controlled flexible ac tra nsmiss 
ion system devices is considered in the paper for the 
first time. In detail, individual generator dynamics are 
characterized with a linear time invariant modeling 
and a nonlinear feedback, with modeling uncertainties 
and power disturbances taken into account. Based on 
such modeling, control laws of SVC's are contrived to 
stabilize the individual generators. The control laws 
of SVC's are developed by exploiting the gain- 
scheduling framework based on the small-gain 
theorem. It is also shown that SVC's can be employed 
for power angle specification in individual generator 
systems. 




FIG. 8 TIME RESPONSES CURVES WHEN THE IGS MAY BE 
UNSTABLE WITHOUT SVC 

The suggested stabilization approach possesses 
several advantages, such as, it can be implemented in 
various decentralized fashions that employ only local 
feedback; or it can cope with damping windings that 
are neglected in the equal area criterion technique; and 
it results in robust stabilization in individual generator 



systems subject to modeling uncertainties and power 
disturbances stemming from governors and loads; the 
SVC control laws can also adaptive to the specification 
requirements on steady-state power angle. It is worth 
mentioning that SVC control laws may also be 
contrived by means of artificial intelligent techniques 
such as fuzzy control and etc. 
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APPENDIX A: PROPERTIES RELATED TO \\G(;K)\\ 

In the main context, we only consider the averagely 
gain-scheduling SVC algorithm that has nothing to do 

with the numeric computation of G(-, . 

However, the small-gain inequality (16) actually can 
have various implementation algorithms beside the 
suggested average one, some of which may need the 

computation of G(-, K)\\ . Hence it is necessary to 
learn more about \\G0, K)\\ to facilitate its numeric 

II 1 1 00 

computation. 

Since ||G(-,X")|| are determined by G(-,X\) , we 
collect basic facts about G(.,J5T ) for a fixed i . Again 

II llco 

noticing that G(-,^ ; ) stands for an LTI system, 

whose transfer function is nothing but G(5,^T ( ) that 
is stable for any K i > , we obtain 

pi;K t )\ = sup \6(j<p,KA 

= sup {(K^Mip 2 ) 1 +(MDIJ) 1 (p 1 Y in ^ 23 ' 

9>e(-oo,oo) 

= { inf FfaK,)}-" 2 



where F(<p,K i ) = (K i -Mp 1 ) 1 + (MD IJ) 2 (p 2 

Taking the derivative on F{jcp,K j ) with respect to (p , 
we obtain 

dF(<p, K. )/d<p = -AM(p(K i - Mcp 1 ) + 2(MD/ jf<p 
which indicates that lower bounds of F{j<p,K j ) may 
appear at 

[>,=o, 

0<K, <M(D/J) 2 /2 
q> x = Q, 9x3 = ±^K t IM -(Dl J) 2 12, 
K. >M(D/J) 2 /2 



284 



International Journal of Automation and Power Engineering (IJAPE) Volume 2 Issue 5, July 2013 



www.ijape.org 



Based on (23) and (24), we can conclude: 

•If 0<^, <M(D//) 2 /2, inf F(<p,K.) is attained 
only at (p x so that 

G(-,/Q =l/K t (25) 

GO 

• If K t >M{DI Jf 12, inf F{q>,K i ) is attained at 
(p x , (p 2 and g) 3 so that 

G(;K,) 



■ max{ 



G(m,^,),go> 2 ,/:,),G(m,^,) 



We have already known that 
Furthermore, we observe 



G(M.*,-) 



=1/^ 



= - M [2K j -M(D/J) 2 ]/(2M)) 2 

+ {MD IJ) 1 [2K I -M{DI J) 2 ]I{2M)Y' 12 

= {{K t -[K l -M(DIJ) 2 /2]) 2 

+ M(D/7) 2 [A: i -M(DIJ) 2 /2])~" 2 

= {(M(D//) 2 /2) 2 +M(D/J) 2 K i 

-M 2 (D/J)"/2y" 2 

= {M (D / J) 2 K i - M 2 (D / J)* /4T" 2 > l/K t 

The last inequality follows from K i > M (D I jf 1 2 2 . 
In summary, it is concluded that the larger K t > is, 



the smaller 



is, and 



G(;K,) 

CO 

0< K, <M(DI J) 2 12, 

II lloo 

K l >M(DIJ) 2 12, 

[M{DIJ) 2 K t -M 2 {DI J) 4 /4}-" 2 >1/ K i 
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